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Abstract 

An accelerated particle sees the Minkowski vacuum as thermally excited, and the particle 
moves stochastically due to an interaction with the thermal bath. This interaction fluctuates 
the particle's transverse momenta like the Brownian motion in a heat bath. Because of this 
fluctuating motion, it has been discussed that the accelerated charged particle emits extra radi- 
ation (the Unruh radiation pQ) in addition to the classical Larmor radiation, and experiments 
are under planning to detect such radiation by using ultrahigh intensity lasers constructed in 
near future [2} [3]. There are, however, counterarguments that the radiation is canceled by an 
interference effect between the vacuum fluctuation and the fluctuating motion. In fact, in the 
case of an internal detector where the Heisenberg equation of motion can be solved exactly, 
there is no additional radiation after the thermalization is completed [HE]. In this paper, we 
revisit the issue in the case of an accelerated charged particle in the scalar-field analog of QED. 
We prove the equipartition theorem of transverse momenta by investigating a stochastic motion 
of the particle, and show that the Unruh radiation is partially canceled by an interference effect. 
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1 Introduction 



Quantum field theories in the space-time with horizons exhibit interesting thermodynamic 
behavior. The most prominent phenomenon is the Hawking radiation [6] and the fundamental 
laws of thermodynamics hold in the black hole background [7]. This indicates an underlying 
microscopic description of the space-time and there are varieties of proposals including those 
based on D-branes [8] or quantum spin foam [9]. For a black hole with surface gravity k at the 
horizon, the temperature of the Hawking radiation is given by 

fin „ H ( M & 



where we have used k — 1/ 4M for the Schwarzschild black hole with mass M . The temperature 
is too small to be observed for astrophysical black holes. A similar phenomenon occurs for a 
uniformly accelerated observer in the Minkowski vacuum [TQl H2J- The equivalence principle 
of the general relativity relates acceleration with gravity. If a particle is uniformly accelerated 
with an acceleration a, there appears a causal horizon, the Rindler horizon, and no information 
can be transmitted from the other side of the horizon. Because of the existence of the Rindler 
horizon, the accelerated observer sees the Minkowski vacuum as thermally excited with the 
Unruh temperature 

= =4x 10 - 23 / [Kl. (1.2) 

2TTck B Vlcm/sV 



Furthermore it is discussed [13] that, if we assign entropy to the Rindler horizon and assume 
thermodynamic relations, the Einstein equation can be derived. 

The Unruh temperature is very small for low acceleration, but the recent development of 
ultra-high intensity lasers makes the Unruh effect experimentally accessible. In the electro- 
magnetic field of a laser with intensity /[W/cm 2 ], an electron can be accelerated to 



a = 2xl0 12 J 2 [cm/s 2 ] (1.3) 

1 W/cm 



and the Unruh temperature is given by 



Tu = 8 x 10~ n W ——— 2 [K]- (1.4) 
y lW/cm 

The ELI (Extreme Light Infrastructure) project [3] recently approved is planning to construct 
Peta Watt lasers with an intensity as high as 5 x 10 26 [W/cm 2 ]. Then the expected Unruh 
temperature becomes more than 10 3 K which is much higher than the room temperature. Can 



1 



we experimentally observe such high Unruh temperature of an accelerated electron in the laser 
field? This is an interesting issue and worth being investigated. 

One such proposal was given by P.Chen and T.Tajima p]. Their basic idea is the following. 
An electron is accelerated in the oscillating electro-magnetic field of lasers. It is not a uniform 
acceleration, but they approximated the electron's motion around the turning points by a 
uniform acceleration. Since the electron feels the vacuum as thermally excited with the Unruh 
temperature, the motion of the electron will be thermalized and fluctuate around the classical 
trajectory (Fig. [T]). Because of this fluctuating motion of an electron, they conjectured that 




X 



Figure 1: Stochastic trajectories induced by quantum field fluctuations. 



additional radiation, apart from the classical Larmor radiation, will emanate. Using an intuitive 
argument, they estimated the additional radiation and called it Unruh radiation. Though the 
estimated amount of radiation is much smaller than the classical one by 10~ 5 , they argued that 
the angular dependence is different. Especially there is a blind spot for the Larmor radiation in 
the direction along the acceleration while the Unruh radiation is expected to be radiated more 
spherically. Hence they proposed to detect the additional radiation in this direction. 

The above heuristic argument sounds correct, but it has been known in a simpler situation 
that such radiation is canceled by an interference effect between the radiation field emanated 
from the fluctuating motion and the vacuum fluctuation of the radiation field [U [5]. The 
cancellation was shown to occuJIj for an internal detector in 1+1 dimensions^). In order to see 



§ We note that, though the radiation vanishes in generic points, a singular behavior of the radiation is pointed 
out to exist on the past horizon [14] . 

' We call a detector with internal degrees of freedom an internal detector in order to distinguish it from an 
accelerated charged particle. In the case of a charged particle, the position of the particle reacts to the thermal 
effect of acceleration. On the contrary, in the case of an internal detector, only its internal degree of freedom is 
excited by the thermal effect. 
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the importance of the interference effect, we briefly sketch the calculation of radiation from the 
internal detector. In these papers [HE], the authors analyzed a uniformly accelerated internal 
detector Q coupled with a scalar field in 1+1 dimensions. The action is given by 

S = S(Q) + 5(0) + e J dr^izir)) (1.5) 

where z{t) = (t(r),x(r)) represents a classical trajectory of the detector. S(Q) and S(</>) are 
quadratic actions of the internal degrees of freedom of the detector (i.e. a harmonic oscillator) 
and the scalar field in 1+1 dimensions respectively; 

S(Q) = jdr Qq(t) 2 - ^ Q 2 (r)) (1.6) 

5(0) = J d 2 x±(d<P(x)) 2 . (1.7) 

Since the coupling term is linear both in Q and 0, the Heisenberg equations of motion can 
be exactly solved. A classical solution of the scalar field is written as a sum of the vacuum 
fluctuation 0^(x), which is a solution to the homogeneous equation in the absence of Q, and 
an inhomogeneous term 4>inh{ x ) as 

0(x) = <j> h (x) + (pinh(x). (1.8) 

The inhomogeneous part 0^ is given by 



<t>inh{x) 



J drG R (x,z(r))^- (1.9) 



where Gr is the retarded Green function of the scalar field. The equation of motion of Q 
becomes 

Q + UQQ = -e— e— — . (1-10) 

dr dr 

Since 0^ is solved linearly in Q(t)., the second term of the r.h.s. of f ll.lOp is linear in Q. It 
can be shown that it becomes a dissipative term 7Q where 7 = e/27r. Hence Q{r) can be solved 
in terms of the homogeneous solution 0/j as 

Q(u) = eh(u)ip(u) (1.11) 

where Q(oS) and (p(oj) are the Fourier modes of Q(t) and 0/ l (z(r)) with respect to r, and 
h{ui) = iu/(uj 2 — Uq — icjj). By inserting (11.111) to (11.91) . the inhomogeneous solution (p inh is 
solved in terms of the vacuum fluctuation 4>h{x). 

Then it is straightforward to calculate the energy-momentum tensor. Since the energy flux 
is written in terms of the 2-point function, they first calculated the 2-point function 

G(x,x') = (0(x)0(O). (1.12) 
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It is written as a sum of the following terms, 

G(x,x') - G (x,x') = ((f>inh{x)<f>inh{x')) + (<Pinh(x)<Ph(x')) + ((f> h {x)<j>inh{x')} (1-13) 

where vacuum fluctuation Gq(x, x') = (4>h(x)(j>h(x')) is subtracted. The first term ((j)inh(x)4 > inh{.x')) 
can be considered as an analog of the Unruh radiation proposed in [I]. It is nonzero because 
the detector is thermally excited from the classical ground state Q — 0. However, Sciama 
et.al. pi] and Hu. et.al. [5] have shown that in (l+l)-dimensional case, the contributions from 
the interference terms (0^(^)0/1 (a 7 )) + (0/i( x )0m/i(a; / )) cancel the radiation (^(1)^/1(2/)) , 
except for the polarization cloud. The polarization cloud is the cloud of radiation field localized 
near the accelerated particle, and does not contribute to the energy-momentum tensor. Since 
we are interested in the energy flux far from the particle in this paper, we do not consider it 
in the following. Hu and Lin [15] extended the calculation to a detector with internal degrees 
of freedom in 3+1 dimensions. The calculation in the 3+1 dimensional case is reviewed in 
Appendix A. 

The physical reason of the cancellation is discussed in [16] and [17] (see also [IS])- The key 
point of their arguments is the invariance of the Minkowski vacuum state under boosts with 
the conservation of energy-momentum tensor. The metric of the Minkowski space is invariant 
under the translation of r. So, after the thermalization occurs, the system becomes stationary 
and the energy-momentum tensor will not depend on r explicitly. In [IT], the following three 
conditions are assumed: 

• Poincare invariance; 

• A detector follows a uniformly accelerated trajectory in flat space; 

• A r-independent coupling. 

With these conditions, the energy-momentum tensor does not depend explicitly on r. They 
have argued that, if the energy conservation holds for the total system, there is no outgoing 
flux. Namely the cancellation of the radiation {4>i n h{x)4>inh{x')) by the interference terms is 
not accidental but a consequence of the boost symmetry. In the case of an accelerated charged 
particle, however, the second condition does not hold since the particle's trajectory is a dynam- 
ical variable and fluctuates. Then the boost invariance of the trajectory is lost. Furthermore, 
the charged particle is accelerated by injecting energy from outside. It is therefore not obvious 
whether extra radiation vanishes or not. 

The purpose of the paper is to investigate a stochastic motion of a uniformly accelerated 
charged particle and to study whether there is additional radiation, Unruh radiation associated 
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with the stochastic motion of the particle. The situation becomes much more complicated than 
the internal detector case because the equations of motion are highly nonlinear. When the 
particle's motion z(t) is affected by the vacuum fluctuation, the Green function Gr(x,z(t)) 
is also changed accordingly unlike the internal detector case. Hence in order to calculate the 
radiation from an accelerated charged particle, we need to approximate the fluctuating motion 
near the classical trajectory and assume that fluctuations are small. We first study a stochastic 
equation for a charged particle coupled with the scalar field [19] . This gives a simplified model 
of the real QED. A self- interaction with the scalar field created by the particle itself gives a 
backreaction to the particle's motion, and it gives a radiation damping term of the Abraham- 
Lorentz-Dirac equation [20]. If we further regard the vacuum fluctuation as stochastic noise, 
the particle's motion obeys a generalized Langevin equation. By solving the Langevin equation, 
we can obtain stochastic fluctuations of the particle's momenta. In this paper, we mainly focus 
on fluctuations in the transverse directions and leave analysis of longitudinal fluctuations for 
future investigation. Some comments are given in Discussions and in Appendix [Bj 

The organization of the paper is the following. In Sec. II, we summarize the basic framework 
of our system, and obtain a generalized Langevin equation for a charged particle coupled with 
a scalar field. In Sec. Ill, we consider small fluctuations in the transverse directions. Then 
the stochastic equation can be solved and we can prove the equipartition theorem for the 
transverse momenta, i.e., a stochastic average of a square of momentum fluctuations in the 
transverse directions is shown to be proportional to the Unruh temperature. We also discuss 
the relaxation time of the thermalization process. Section IV is the main part of the paper. 
We calculate radiation emitted by a charged particle in the scalar-field analog of QED. The 
interference terms partially cancel the radiation coming from the contribution {4>i n h{x)4>inh( x ')) > 
but unlike the case of an internal detector, they do not cancel exactly. In Sec. V, we obtain 
a similar stochastic equation for an accelerated charged particle in the real QED, and show 
the equipartition theorem for transverse momenta. Section VI is devoted to conclusions and 
discussions. In Appendix A, we review the calculation of radiation for a uniformly accelerated 
internal detector in (3+1) dimensions [15]. In Appendix (Bj we consider fluctuations in the 
longitudinal and temporal directions. We use Planck units c = H — ks — 1 in most of the 
paper. 

2 Stochastic equation of an accelerated particle 

We consider the scalar-field analog of QED. The model is analyzed in [19] and here we briefly 
review the settings and the derivation of the stochastic Abraham-Lorentz-Dirac (ALD) equa- 
tions. In |19| . the authors used the influence functional approach, but here we take a simplified 
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method. The system composes of a relativistic particle z M (r) and the scalar field <f)(x). The 
action is given by 

S[z, 0, h\ = S[z, h] + S[(j)] + S[z, 0], (2.1) 

with 



S[z, h] = 


—m 


S[*] = 




S[z,</>] = 





4. 1 



(2.2) 

- gv-,.^-,, , (2-3) 

d*x j(x;z)(f)(x). (2.4) 
The scalar current j(x; z) is defined as 

j(x; z) = e I dr y/Ji%5 4 (x - z(r)), (2.5) 



where e is negative for an electron. We can parametrize the particle's path satisfying z 2 = 1 
by taking r properly. 

The equation of motion of the particle is given by 

= F»- [ d A x 53 } X ) Z } <p(x) (2.6) 
J oz^(t) 

where we have added the external force so as to accelerate the particle uniformly: 

= ma(i 1 , i°, 0, 0). (2.7) 

Then a classical solution of the particle (in the absence of the coupling to <fi) is given by 

1 l 

z^ = (- sinh ar, - cosh ar, 0, 0). (2.8) 
a a 

Note that the external force satisfies F^z^ = and therefore the classical equation of motion 
preserves the gauge condition z 2 = 1. From the definition of the current (I2.5p . it is easy to 
prove the identity, 

/ d ' X l^y f{Xj = e ^/(*)U M (2-9) 

where at M is given by 

at M = z u z [v d A - z^. (2.10) 

Here we have used the gauge condition z 2 = 1 and z ■ z — 0. Hence the equation of motion 
( 12. 61) becomes 

mi* = F" - efi^zlj)) (2.11) 
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Since the differential operator 73* M satisfies i M at M = for a classical path satisfying the gauge 
condition, the stochastic Eq. (12. lip continues to preserve the condition z 2 — 1. The second 
term of (12.111) represents a self-interaction of the particle with the radiation emitted by the 
particle itself. 

The equation of motion of the radiation field d^d^ix) = j(x) is solved by using the retarded 
Green function Gr as 

00) = (j) h (x) + <f) inh , 4> inh = / d 4 x'G R (x,x')j(x'; z) (2.12) 



where <ph is the homogeneous solution of the equation of motion and represents the vacuum 
fluctuation. The retarded Green function satisfies 

d^d^Gnix, x') = 5 {4) (x - x') (2. 13) 

and is given by 

G R {x,x') = i(l<P(x)^(x')))9(t-t') (2.14) 
= 6(t - t')5((x - x'f) = m-t')-r) 

where r 2 = |x — x'| 2 . Inserting the solution f |2.12j) into (I2.1ip . we have the following stochastic 
equation for the particle 



V) = F»(z(t)) - U h (z(r)) + e J dr' G R (z(r), z(/))\ . (2.16) 



Here we have used the gauge condition z 2 = 1. The operator at M acts on z(r). The homogeneous 
part (f>h(z(T)) of the scalar field describes Gaussian fluctuations of the vacuum, and hence the 
first term in the parenthesis can be interpreted as random noise to the particle's motion. 
Expanding (fi h as 

^-IS-^ a ^'' +a ^ (2 - i7 » 

the vacuum fluctuation is given by 

= -^ (t _ t ,_\ l€ y_ r r (2.18) 

It is essentially quantum mechanical, but if it is evaluated on a world line of a uniformly 
accelerated particle x = z(t),x' = z(r'), it behaves as the ordinary finite temperature noise. 

The second term in the parenthesis of (12.161) is a functional of the total history of the 
particle's motion z(t') for r' < r, but it can be reduced to the so called radiation damping 
term of a charged particle coupled with a radiation field. It is generally nonlocal, but since 
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the Green function damps rapidly as a function of the distance r, the term is approximated by 
local derivative terms. First define ?/ M (s) = z M (r) — z^{t') where r' = r — s with r kept fixed. 
Then it can be expanded as 



s 2 .... , , s 3 



y»(s) = sz^r) - -F(r) + -z"(t) + ■ ■ • . (2.19) 
2 o 



A square of the space-time distance a is given by 



cr(s)=y% l = s 2 ^l-^(z) 2 + --- ) (2.20) 
and 

^ = 2y»y, = 2s (l - S j(-zf + -•). (2.21) 

In deriving them, we have used the gauge condition (i) 2 = 1, z • z = and z • 'z = — {z) 2 - The 

derivative <9 M appearing in the operator ut M can be written in terms of 4-, when it acts on a 
function of a, as 

^ = ¥-4- = 2y M f^V 1 # = (2.22) 

0za« da y ^\dsJ ds y v y v ds v ; 

z, - S -z, + S -(z, + i^z) 2 ) + ■ • • ) ^- (2.23) 
Hence the second term in the parenthesis of (12. 161) can be simplified as 

f dr'G R (z(r),z(T f )) (2.24) 

J — oo 
/>oo 

; 2 / ds^ fl (s) (2.25) 
Jo 

jf U(t)G r (s) + a^^^G?^) + (i^ 2 + z ,)- ^G R ( S ) + ( S 3)j . (2.26) 



e 2 " 



--e 2 



In the first equality, we have neglected the singular term proportional to 6(a). The first two 
terms can be absorbed by mass renormalization. The s integrals of them are divergent, so we 
assume that the renormalized mass becomes finite after the mass renormalization. The last one 
is the radiation reaction term and can be evaluated by using the identity 



dss 2 —G H (s) = dss 2 —^ = (2.27) 

ds K 1 J dsAns 2vr v 1 

After the mass renormalization, we get the following generalized Langevin equation for the 
charged particle, 

p 2 

m? — (z^z 2 + "if") = -eat"^*). (2.28) 

127T 



This is an analog of the Abraham-Lorentz-Dirac equation for a charged particle interacting with 
the electromagnetic field. The dissipation term is induced through an effect of backreaction 
of the particle's radiation to the particle's motion. Note that, if the noise term, the r.h.s., is 
absent, the classical solution (12.81) is still a solution to Eq. (12. 28 p . 

3 Thermalization of transverse momentum fluctuations 

The stochastic Eq. (I2.28P is nonlinear and difficult to solve. Here we consider small fluctuations 
around the classical trajectory induced by the vacuum fluctuation <ph- Especially, we consider 
fluctuations in the transverse directions perpendicular to the direction of the acceleration. The 
classical trajectory of the particle does not move to the transverse directions, and we can easily 
treat fluctuations in these directions. In contrast, the particle is accelerated strongly to the 
longitudinal direction and it is more difficult to separate fluctuations from the classical solution. 
We shortly discuss longitudinal fluctuations in Appendix [0 

First we expand the particle's motion around the classical trajectory Zq as 

z'>(t)=zS + 6z> 1 . (3.1) 

The particle is accelerated along the x direction. Now we consider small fluctuation in transverse 
directions z l . By expanding the stochastic Eq. ( 12. 28ft . we can obtain a linearized stochastic 
equation for the transverse velocity fluctuation 8v l = 5z l as, 

e 2 

mSv l = ed^ h + (5v l - a 2 5v l ). (3.2) 

127T 

Performing Fourier transformation with respect to the trajectory's parameter r 

Sv\r) = J ^-5v\u)e-^, ^ fc (r) = J ^d lV (co)e~^, (3.3) 
the stochastic equation can be solved as 

5v\u) = eh(uj)diip(uj) (3.4) 

where 

h ( w ) = — — , e^wr ( 3 - 5 ) 

—imUJ H ^ 

The vacuum 2-point function on the classical trajectory can be evaluated from ( 12.181) as 

(diM^djMx'Vl^TU^ziT') = ^ (( t _ t /_^ £ ) 2 _ r 2) 2 ( 3 - 6 ) 

=_ °J1 (3 7) 

32?r2 sinh 4 ( "(t-t'-^ ' 
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It has originated from the quantum fluctuations of the vacuum, but it can be interpreted as 
finite temperature noise if it is evaluated on the accelerated particle's trajectory [10] 0. By 
making the Fourier transformation with respect to r, we have 

{d i tp{u)djtp{u')) = 2tx5(u + (3.8) 

where 

- 32tt 2 7_„ ^(2^) "6tt 1 - e -W«" ^ 
For small u, this is expanded as 

I( u ) = (a 2 + a-KUJ + ■■■). (3.10) 



127T 2 

The expansion corresponds to the derivative expansion 

a 2 

(<9^(x)<9^ h (£0)U= 2 (T) i:C '=^T') = -^^(r ~ r ') ~ i 7o~V( T ~ r ') H " ( 3-11 ) 

If we approximate the 2-point function by the first term, the noise correlation becomes white 
noise. The coefficient determines the strength of the noise. We show that it is consistent with 
the fluctuation-dissipation theorem at the Unruh temperature. 

By symmetrizing it, i.e., (d<f)(x)d<j)(x'))s = {{d<fr(x),d<j)(x')})/2, the correlation function 
becomes 

W = — cothl — ). (3.12) 

It is an even function of uj. The correlators I(oS) and Is{w) should be regularized for large u 
where quantum field theoretic effects of electrons become important. 

The expectation value of a square of the velocity fluctuations in the transverse directions 
can be evaluated as 

(6v\t)Sv>(t')) s = e 2 j d ^{d i y{u)d i y{u')) s h{u)h{u')e^ + ^ (3.13) 
= e% [ ^l s {u)\h{u)\ 2 e-^-^ (3.14) 



^ e 2 5 .. / _^ °L e -Mr-r') (315) 

24vr 3 {mujf + (^) {u 2 + a 2 ) 2 V ' ' 



"The finite temperature (Unruh) effect is caused by the appearance of the horizon for a uniformly accelerated 
observer in the Minkowski space-time and analogous to the Hawking radiation of the black hole, but you should 
not confuse the radiation we are discussing in this paper with the Hawking radiation. The accelerated observer 
sees the Minkowski vacuum as thermally excited, but it is excited from the Rindler vacuum (not from the 
Minkowski vacuum) and the energy-momentum tensor remains zero as ever. The radiation discussed in the 
paper is, if exists, produced by an interaction with the vacuum and the accelerated charged particle. 
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The denominator has four poles at u = ±zf2± where 

n + = ^(l + 0(o a /m 2 )), (3.16) 

fi_ = ^(l + 0(a 2 /m 2 )). (3.17) 

The acceleration of an electron in high-intensity laser fields in the near future can be at most 
0.1 eV and much smaller than the electron mass 0.5 MeV. Hence, the values of these poles 
satisfy the following inequalities, 

fi+>a>fi_. (3.18) 

Since the energy scale of the dynamics of the accelerated particle is much smaller than the 
electron mass, the poles at ±iQ + should be considered spurious and we should not take the 
contributions of the residues at ± «ii By taking the residue at ±oj = we can evaluate 
the integral and get the following result, 

777 1 fl H 

-(5v\r)5v^r)) = -—5 %3 (l + 0(a 2 /m 2 )) . (3.19) 

Here we have recovered c and h. This gives the equipartition relation for the transverse mo- 
mentum fluctuations at the Unruh temperature Tjj = ah/2irc. The typical energy scale of 
the fluctuating motion is given by the value of the pole Since it is much smaller than 
the acceleration, the derivative expansion with respect to u/a is justified for the transverse 
fluctuations. 

The thermalization process of the stochastic Eq. ( 13. 21) can be also discussed. For simplic- 
ity, we approximate the stochastic equation by dropping the second derivative term. This 
corresponds to the derivative expansion with respect to u/a. Then it is solved as 

SvHt) = e~ n - T 5v l (0) + — [ T dr' 9 4 0(^(r , ))e- Q - (r " T ' ) . (3.20) 

m Jo 

The relaxation time is given by tr = l/f2_. The momentum square can be calculated as 

(^(r)^(r)) =e~ 2n - T 5v\0)5v j (0) 



2 

e 



^0 



dr' / dr" e- ^7 - (r - T ' ) e" ^ - (T " r " ) (^(^(r / ))a i 0(z(r /, ))) 



- 2Q -W(0)^'(0) + - e~ 2n - T ). (3.21) 



** Or we can simply approximate the denominator by dropping the u term. Then only the poles at ±ifL 
survive. 



11 



For r — > oo, it approaches the thermalized average (13.191) . The relaxation time in the proper 
time can be estimated, for a = 0.1 eV and m = 0.5 MeV, to be 

12-Kirt r 
TR = —— = 1.4 x 10" 5 sec. 3.22 
a z e z 

This relaxation time should be compared with the laser frequency. The oscillation period of the 
laser field at ELI is about 3 x 10~ 15 seconds and much shorter than the above relaxation time. 
Hence an accelerated electron by the ELI laser does not thermalize during each oscillation and 
the uniform acceleration is not a good approximation. Even in such a situation, if an electron is 
accelerated in the laser field for a long time, it may feel averaged temperature. In order to fully 
understand the Unruh effect of nonuniform acceleration, we need to investigate the transient 
phenomena. It is beyond the present calculation, and we leave it for future investigations. See 
e.g. [H] for the Unruh effect of the circular motion, . 

The position of the particle in the transverse directions fluctuates like the ordinary Brownian 
motion in a heat bath. A mean square of the fluctuations is given by 

R\r) = ((Ar) - ^(0)) 2 ) = 2DU- 2Q + ) ■ (3.23) 

i=y,z 

The diffusion constant D is given by 

= ^ = 4. (3-24) 
iZ_m ae z 

which is estimated for the above parameters as D ~ 7.8 x 10 4 m 2 /s. In the Ballistic region 
where r < tr, the mean square becomes 

fff T ) = (3.25) 
m 

while in the diffusive region (r > tr), it is proportional to the proper time as 

R 2 (r) = 2Dt. (3.26) 

As the ordinary Brownian motion, the mean square of the particle's transverse position grows 
linearly with time. If it becomes possible to accelerate the particle for a sufficiently long period, 
it may be possible to detect such a Brownian motion in future laser experiments. 

4 Quantum Radiation by Transverse Fluctuation 

Once we obtain the stochastic motion of the accelerated particle, it is straightforward to calcu- 
late the energy flux of the radiation field emitted by this particle. In this section, we calculate 
the radiation induced by the fluctuation in the transverse directions. 
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4.1 2-point function 

First we evaluate the 2-point function 



G{X,X') - G Q {X,X') =((t)inh(x)<fiinh(x')} + (4>inh{x)4>h{x')) + {(j) h {x)(f)i nh (x')) . (4.1) 

The inhomogeneous solution <p in h is produced by the accelerated charged particle while the 
homogeneous solution 0^ represents the vacuum fluctuation of the quantum field 0. The Unruh 
radiation estimated in [TJ corresponds to calculating the 2-point correlation function of the 
inhomogeneous terms (4>i n h{x)4>i n h{x')) ■ As we will see later, this term also contains the classical 
Larmor radiation. However, this is not the end of the story. As it has been discussed in [I], the 
interference terms {4>i n h4>h) + {(frh&nh) cannot be neglected and may possibly cancel the Unruh 
radiation in (4>inh4>inh) after the thermalization occurs. This is shown for an internal detector 
in (1+1) dimensions, but it is not obvious whether the same cancellation occurs for the case of 
a charged particle we are considering. 

The inhomogeneous solution of the scalar field is written in terms of the position of the 
accelerated particle z^(t) as 

(pinh(x) =e J drG R (x - z(r)) (4.2) 

I >( t - AT)W (,-^))') = _^ 

J 2tt 4vrp v ; 

where p is defined by 

P=z{t x 1) ■ (x-z(t*)). (4.4) 

Because of the step and the delta functions in the integrand of (14.31) . r'f_ satisfies 

(x - z(t x 1)) 2 =0, x° > z°(rf ), (4.5) 

which is the proper time of the particle whose radiation travels to the space-time point x. Hence 
z(tH) lies on an intersection between the particle's world line and the light-cone extending from 
the observer's position x (see Fig. 2). We write the superscript x to make the x dependence of 
r explicitly. The meaning of the subscript (— ) will be made clear later. By using the light-cone 
condition, p can be rewritten as 

P(x) = ri - ^ (4.6) 



dr 

where v = r(rf) = x — z(r^) and r = |r|. It is the spacial distance for the observer moving 
with the particle. 
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The particle's trajectory is fluctuating around the classical trajectory and can be expressed 
as z = zq + 5z + 5 2 z + ■ ■ ■ where the expansion is given with respect to e. In terms of the 
expansion of the spacial distance p = p + 5p + 5 2 p + • • • , the inhomogeneous solution (14. 3 j) 
becomes 



■'inh 



l- S ± + ( S ±\ -*£ + ...!. (4.7) 



4vrp \ Po VPo/ Po 



The first term is the classical potential, but since the particle's trajectory deviates from the 
classical one, the potential also receives corrections. Here po, Sp and S 2 p are given by 

po =i (r*) ■ (x - z (t*)) (4.8) 
5p =5z(t x _) ■ (x - z (r x )) - z (r x _) ■ 5z(t x _), (4.9) 
S 2 p =S 2 z(t x ) ■ (x - z (t x )) - 5z{t x _) ■ Sz(t x ) - z (r x ) ■ 5 2 z(t x ). (4.10) 

From now on we only consider the transverse fluctuations. Then z ■ 5z = is satisfied. As 
seen from (14.61) . p is proportional to the spacial distance from the particle to the observer. The 
variation of p becomes negligible for large distance r if we take a variation of (x — Zq{t)) in 
p. On the contrary, if we take a variation of zq, Sp or 5 2 p is still proportional to the spacial 
distance r. Hence for large r, we can approximate the variations by 

5 P ~6z(tZ)-(x-z (t*)), (4.11) 
5 2 p~5 2 z{t x ) ■ {x-z {t x )). (4.12) 

Note also (5 2 z l ) = since the velocity in the transverse directions fluctuates uniformly and its 
expectation value vanishes. 

4.2 Correlations of the inhomogeneous terms 

Now we calculate the 2-point function explicitly. If we take the classical part without the 
fluctuation of p, the 2-point function becomes 

f e \ 2 1 

G(x, x) - G (x, x) (<j) in h(x)<j>i n h(x')) = — —rr — -r-^- (4.13) 

V4vr/ p (x)po{x') 

This gives the classical radiation corresponding to the Larmor radiation. The interference term 
vanishes because the 1-point function vanishes identically (<ph) = 0. 

Corrections to the above classical Larmor radiation are induced by the transverse fluctuating 
motion 5p. First we consider the 2-point correlation function between the inhomogeneous terms 
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up to the second order of the transverse fluctuations. Since (S 2 p) = 0, we have 

'e_V 1 ( 1 + (Sp(x)5p(y)) + ((W) 2 ) + ((Ml/)) 2 ) 
An) p (x)p (y) V p (x)p (y) Po0*0 Poiv) 



(4.14) 



Note that all the terms in the parenthesis behave constantly as the distance r between the 
observer and the particle becomes large. The first term gives the Larmor radiation mentioned 
above. The other terms correspond to the radiation induced by the fluctuations. Its calculation 
is easy, because one can write (Sp5p) in terms of (8z l 8z l ) = (5v l 5v t ) which we have already 
evaluated in the previous section. With the expression (13 .151) . it becomes 



(4>inh{x)4>inh{y)) 



e y 1 
4tt/ p (x)p Q (y) 



e 2 / ^\h{u)\H{u) 



x 



x % y l e luJ ( T - r -) x l x % y l y l 



po(x)p (y) po(x)p (x) po(y)po(y) 



(4.15) 



As before, since we are considering the fluctuating motion whose frequency is smaller than the 
acceleration, we may as well approximate I(uj) by a 3 /127r 2 . In order to calculate the symmetrized 
correlation function between x and y, I(cu) is replaced by Is(oj). 

4.3 Interference terms 

Next let us calculate the interference terms. They are rewritten as 

{<p inh {x)<p h {y)) + {<p h {x)<p inh {y)) = -— ^ + j . (4.16) 

Calculation of the interference terms are more complicated since we need to evaluate the fol- 
lowing correlation functions: 

(dp(x)My)} = - x*(5z\T*)My)} (4.17) 

= -ex i f pe"^-h(u)(d tV (u)My)) (4-18) 
J ^ 

(<f> h {x)8p{y)) = - y'iMWz'iil)) (4.19) 

= -ey l j ^e^-h{u){<t> K {x)di<p{u)). (4.20) 

Since two terms (di<p(oj)<j)h(y)) and (<f>f l (x)di<p(u)) are related by 

(dM^Mv)) =({Mv)dM-u>W, (4.21) 
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it is sufficient to calculate one of them. From the definition of ip in (13.31) . the interference term 
{(j)h{x)dip{uj)) is written as 

(Mx)dMu>)) = j dre^ (^-(^{x)^))^ ^ ^ (4.22) 

( „ .,' ,-> ->,,, ) (4.23) 



y=z(r) 

^ r 7Zn ..n/_\ 77T5 TZi I77Z\v3> ZiT> (4.24) 



4tt 2 7 V ^ (x° - y° - ze) 2 - - ^) 2 

1 g 

"i^ 2 ^ 7 (^-z^r)-^) 2 -^ 1 -^)) 2 -^' 

where = (x 2 ) 2 + (x 3 ) 2 is the transverse distance. We first evaluate the integral and then 
take a derivative. The integral 

/ e iuiT 
dr (xO-,o (r) _, £)2 _ (x i_,i (r))2 _ xr ( 4 - 25 ) 

can be evaluated by the contour integral in the complex r plane. The positions of the poles are 
given by a series of points 

m 2mri 

r± =T± + ie, (4.26) 

(X 

where n is an integer. T± are complex numbers whose imaginary parts are or ir/a and satisfy 



a o /_ , 4 



Here we have defined 



" 2 - - x* 

„0 „1 „. _0 1 1 



L — — x% + — , (4.28) 

or 



u=x u -x\ v = x {i + x\ (4.29) 

Note the relation e aT+ e aT ~ = —v/u. The positions of the poles reflect the finite temperature 
property of the uniformly accelerated observer. In the following we will consider two different 
types of observers as shown in Figj2j The first observer is to observe the radiation in the 
right wedge (Or) while the second one is in the future wedge (Op)- For both cases, v > is 
satisfied and the radiation can travel causally from the particle to the observers. There are 
two different types of poles t±. A pole at r_ = T_, which is real, is located at a classically 
acceptable point. Namely, r_ is the proper time of the particle whose radiation travels to the 
observer in a causal way. The other pole at T + is more subtle. For u < (in the right wedge), 
T + = r+ is real and corresponds to the advanced causal proper time. For u > (in the future 
wedge), T + = r+ + in / a has an imaginary part and one can interpret it as the proper time of 
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Figure 2: The hyperbolic line in the right wedge denotes the world line of the particle. The 
points Of and Or are observers in the future and right wedges, respectively. For an observer 
in the right wedge, the light-cone of the observer has two intersections with the world line, 
and the proper time of the intersections is given by t±. For an observer in the future wedge, 
there is only one intersection on the particle's real trajectory which corresponds to t^. The 
other solution Tj = + in /a is complex. One may interpret this complex proper time as the 
intersection between the light-cone of the observer and the world line of a virtual particle with 
a real proper time r+ in the left wedge. The superscript letters R or F are used to distinguish 
two different observers, but we do not use them in the body of the paper to leave the space for 
the observer's position x. 
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a trajectory of a virtual particle in the left wedge, as in Fig. [2j In the following, we drop the 
superscript F or R. In the region where v < 0, 4>i n h does not exist and no nontrivial correlation 
is observed there. 

The residue of the pole at t± is given by -e^ T ±/(2p(r£)) where 

p{rl) = z{rl) ■ (x - z Q {rl)) (4.30) 
= \(ue aT ± +ve- ar "). (4.31) 

Because of the periodicity, p{r±) is independent of n. The integral is now given by 

— IT 7 

P M =—^- [ (e^--e^lZM), (4-32) 

where 

Zju) = e wui/a e{u) + 9{-u) (4.33) 
Po — p( T ~) can be rewritten in terms of L 2 as 



Po = ^L±+^uv. (4.34) 

Note that the relation p(r?) = — po follows the identity e aT+ e aT ~ = —v/u. 

The second term of the parenthesis in ( 14.321) depends on r + . With naive intuition based on 
classical causality, the term may be removed by hand, but the calculation of the interference 
terms is essentially quantum mechanical, and it should not be neglected. It is puzzling how we 
can physically interpret such r + dependence of the integral. 

Taking a derivation of P(x,u), we obtain ((f>h(x)di<p(uj)) as 

(««)M<-)> =i^3^T((^ + + + -) e + ^^>> (4 ' 35) 



Here we have used the following identities, 

dp a 2 L 2 
dx i 2p ' 



(4.36) 



£f = ± ix*, (4.37) 

where z is the transverse direction. The second identity can be obtained by differentiating 
(x — z{t±)) 2 = with respect to x\ see (14.49 p . 
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The whole interference terms are now given by 



(4> h {x)4>inh{y)) + (<Pinh{x)4> h {y)) 



—%ae 1 x % y % 



(4vr) 2 po(x) 2 p(z/) 2 
+ e- iu ' {T +- T - ) h(-u) 



du 



2?r 1 - . 

2 Po (x) 



-2-rruj/a 

iu 
a 



h(-u) 



all 



ZJ-u) 



2p (x 



1U 

a 



-h(u) 



aL 2 iuj 
+ — 



2po(y) 



a 



IU 



2p (y) 



Zy(-Uj) 



(4.38) 



4.4 Partial cancellation 



In the following, in order to see whether there is a cancellation between the interference terms 
and the correlation function of the inhomogeneous terms, we look more closely at the first 
term in the parenthesis of ( 14.381) which depends only on r_. Note that the correlation function 
of the inhomogeneous terms (14. 15|) depends only on r_, and the r + depending terms in the 
interference terms cannot be canceled with the correlations of the inhomogeneous terms. 
Note that, by using the relation 



h(u) + h(-u) 



[cu 2 + a 2 )\h{u)\< 



one can show that a part of the interference terms in (14.381) 

doj 1 



iae 2 x l y l 



(4ir) 2 p (x) 2 p (y) 2 
can be rewritten as 

e \ 2 x l y l 



2tt 1 



h(-u) h h(oo) — 

a a 



An J po(x) 2 p (y) 2 J 2tt 



du 



^e- i ^ T -- T -)\h(u)\ 2 I{u) = - (-- 



e \ 2 (Sp(x)Sp(y)) 



4tt/ po(x)po(y) 



(4.39) 



(4.40) 



(4.41) 



This cancels the first correction term in the correlations of the inhomogeneous parts in (14.151) . 
Note that this canceled term was obtained by taking a derivative of e tuJT - in P(x,u). 

We have seen that there is a partial cancellation between the interference term and the 
correlations of the inhomogeneous terms, but the other terms are not canceled each other. 
Then, summing up both contributions, (14.151) and (14 .38 p . we get the following result of the 
2-point function; 

„2 



{<j>{x)(j>{y)) - (<j) h (x)(f) h (y)} 



(An) 2 p (x)p (y) 



F(x,y) 



(4.42) 
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where 

T-,t \ ■> 2 f dw \h(oo)\ 2 Tl . / / x l 
F{x, y) = l + e 2 / — I(u) 1 



2tt 6tt \\p (x)J \po(y) 



ia x l y l [duo 1 



Po{x)p (y) J 4nl-e~ 2 ™/ a 



Po{x) po(y) 



\p (x) a* J V Po(y) a 2 / 

(4.43) 

The first term in F(x, y) is the classical effect of radiation corresponding to the Larmor radia- 
tion. The second term comes from the inhomogeneous term ((5 p(x) / po(x)) 2 ) + ((6 p(y) / Po(y)) 2 ) ■ 
The third term comes from the interference term, which is obtained by taking a derivative of 
p(x) in P(x,cj). The fourth term is also an interference effect and depends on r + . 

Let us compare the above result with the calculation for an internal detector. In the case of 
an internal detector in (1+1) dimensions, the radiation is canceled by the interference effect, and 
there are no terms depending on r + . In the case of an internal detector in (3+1) dimensions, 
there are r + -dependent terms. But if we neglect these terms, it was shown [19] that the 
interference terms completely cancel the radiation. The calculation is reviewed in Appendix A. 
In the case of a charged particle, since the position of the particle is fluctuating, only a part 
of the terms is canceled. In the following we focus on the r + -independent terms because the 
rapidly oscillating function e~ lulT remains even after setting x = y and suppresses the u integral 
of the r+-dependent terms. 

For a symmetrized 2-point function, F(x,y) is replaced by F$(x,y) 

. du\h(u)\ 2 T . J f x l \ 2 f y l 
2vr 6vr \\p (x)J \po{y), 

- f -coth (™) e-^--n) ( K-u,)* - h{u)^t\ 

Po(x)po(y)J 8tt \aJ V Po(x) V ' f>o{y) ) 

+ r + -dependent terms. (4.44) 
4.5 Energy- momentum tensor 

In the remainder of this section we consider the radiation emitted by the accelerated particle. 
The energy-momentum tensor of the scalar field is given by 

(Tfiv) = (: d^drf - ^<9 a <R0 :) s . (4.45) 
Hence we can evaluate it by taking a derivative of the 2-point function (14.441) . 
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The following relations are useful in taking derivatives: 



d^po =(i'o • (x - zo) - l)fyr_ + z 0fl (4.46) 
= - d,r. + z , (4.47) 

Po 

In the last line of the first equation, we used the explicit form of the classical solution (12.81) 
and zq ■ x = —a 2 L 2 /2. The derivative c^r_ was obtained by taking a variation of the light-cone 
condition (x — z q (t*)) 2 = 0; 

X 

2{x v - z 0v )(5x v - z»5t x _) = — > = ^i^. (4.49) 

ox po 

In particular, u and w derivatives are given by 

o u r- = — , o v t- = — (4.50) 

zp zpo 

a 2 L 2 av a 2 L 2 1 
duPo = 7T~ d uT- + — -, d v p = ^-d v r- + — — (4.51) 

where u z = — e~ aT ~ /a, v z = e aT ~ /a. From ( I4.48p . we have (<9po) 2 = a 2 x 2 . Since {x — z{r)) 2 = 0, 
x 2 ~ 0(r) and (<9po) 2 is approximately proportional to the spacial distance r, not r 2 . On the 
other hand, since L 2 = — x 2 + 1/a 2 is 0(r), <9 M po itself is growing as 0(r). 

First we calculate the classical part of the energy-momentum tensor. It becomes 

e 2 {d fl p d l/ p - 9 -fd ol p d a p () ) 



e 2 d^p Q d u p Q 



(4.53) 



Note that d a pod a po does not make a contribution here, since it is of the order of po & t infinity 
while d^podvPo is in general of order p\. This part of the energy-momentum tensor corresponds 
to the classical Larmor radiation and behaves as 1/p 2 , ~ \jr 2 at infinity. The term io^(rf) 
in d^po seems to be negligible, since it is 0(1) while <9 M po is 0{r). However, care should be 
taken because ioju( r -) = (cosh ar^, sinh arf , 0, 0) behaves singularly if the observer is near the 
horizon. 

Next we evaluate the other parts of the energy-momentum tensor. We especially consider 
the (u,u) and (v, v )-components in the following. From (14.441) . extra terms of the energy- 
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momentum tensor besides the classical ones are given by 



i\2 



Pi 



6ma 2 IiL 2 



K Po 
2mh 



T, 



e 2 a 2 L 2 



(ml 3 d^d v r a _ 



e 2 I m 



^ L2 ^-»ru ■ ~-A.ru, , 127rL2 



e 2 I 



247T P o^ M 

where we have defined the following u; integrals 

f did., , ,,, , /7TW 

* = y te |ftM| c ° th (t 

/ s =y-iM-)i 2 coth(- 



(4.54) 



J TO = y ^|/i(w)| 2 coth (™) (w 3 + aV) 



= h + a 2 h. 

These integrals can be similarly evaluated as in Sec. Ill, and we have 

/, 3 



I m ~ 



2mae 2 ' 
a 2 h. 



(4.55) 

(4.56) 

(4.57) 
(4.58) 



(4.59) 
(4.60) 



Because of the inequality f2_ <C a, terms containing J 3 are generally negligible compared to 
other terms; I3 ~ f2 2 I\ <C a 2 I\. 

Near the past horizon, the v — > 0, the -u-derivatives of po and become very small and 
negligible. On the other hand, w-derivative of r_ becomes potentially large, w-derivatives of 
them are approximately given by 



av 



d v p -> 



au 
Y' 



(4.61) 



A singular term of d v po near v ~ is canceled and it remains finite near the past horizon. Hence 
the second term in (14.541) proportional to (<9t_) 2 may becomes large there. However, there are 
two reasons that the term cannot grow so large. One is a suppression by the u integral, which 
is proportional to a very small coefficient I3. The other reason is the overall factor (x 1 ) 2 /pq. 
Since the observer is much further than the acceleration scale 1/a from the particle, L 2 is much 
larger than 1/a 2 . Then p = (a/2)^L 4 + (A/a 2 )uv can be approximated by po ~ ( a /2)|^ M | 2 
and (x 1 ) 2 /pq is also suppressed. Because of these two reasons, the singular behavior near the 
past horizon seems to be difficult to be observed experimentally. 
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5 Thermalization in Electromagnetic Field 



In this section, we consider the thermalization of an accelerated charged particle in the elec- 
tromagnetic field. Calculations of the energy-momentum tensors are more involved and left for 
a future investigation. We study the thermalization of the transverse momenta of a uniformly 
accelerated particle in an electromagnetic field. The calculation is almost the same, but due to 
the presence of the polarization, several quantities become twice as large as those in the scalar 
case. 

The action is given by 

Sem = -m J dr^% - J d 4 x f(x)A ll (x) ~\Jd A x F^F^, (5.1) 
where the current is defined as 

f(x) = e J dTz^(r) 5\x - z{r)). (5.2) 

The equations of motion are 

mz^ =eF flu z u (5.3) 
d^(x) =f. (5.4) 

Using the gauge 

WAp = 0, (5.5) 

the equation of motion for A^ becomes 

dPdrA" = f. (5.6) 

One can solve this equation as 

Af, =A hll + J d 4 y G R (x,y)j fJ ,(y) 

=A hfl + e J drG R (x,z(r))z^T), (5.7) 

where A^ is the homogeneous solution of the equation of motion which satisfies d 2 A^ = 0. 
Gr{x — y) is the retarded Green function 

G R (x, y) = 9(x° - y°) * ((g ~ y)2) , d 2 G R (x, y) = 5\x - y). (5.8) 
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Inserting the solution of A^x) back to the equation of motion for z**, we obtain the following 
stochastic equation 

mz^r) =F ll + e(d^A hu (z) - d v A hl]L {z))z v 

+ e 2 f dr'z^r)(z v (r')d, - z^r')d v )G R (z(r) : z(r')). (5.9) 



The second line is the radiation reaction which can be treated similarly to the scalar case. It 
becomes 



e 2 



Jdr' (?z {v d A -j(z%-z^d u )\^-G R (z(T),z(r')) 
e 2 r ds S ^z»{T) + z (1 {T)?{T)) d 6{S2) 



3 v MW " w v " ds 2ix 

~(z^ + z,z 2 ), (5.10) 

07T 

which has exactly the same form as the scalar case, but the coefficient is twice as large since the 
gauge field has two different polarizations. This is the Abraham-Lorentz-Dirac self-radiation 
term. 

For small transverse momentum fluctuations 6v % = 5z l , we can simplify the stochastic 

equation similarly to the scalar case in previous sections. It can be solved in terms of the 
homogeneous solution of the gauge field as 

5v\u) = -eh^iv^di + 5> • d))A a h , (5.11) 

where 

h(uj) = I ■ (5.12) 

—imui + + a ) 

The noise correlation of A^ in the r.h.s. of (15. lip can be evaluated as 

M + 5i(v ■ d))(v' ^ + aj(t/ ■ ff)){A%(z)AlW) = 6 f s — Y . (5.13) 



2 

It is also twice as large as the scalar case. Note that the quantity is gauge invariant 

{z a k„ - 7j a „(i ■ k))%K ~ Vfivtf ■ k'))k a k'P = 0. (5.14) 

Hence performing similar calculations to the scalar case, the fluctuations of the transverse 
momenta become 

"M«))^f,(w(J). (5.15) 

Since the coefficient of the dissipative term is twice as large as the scalar case, the relaxation 
time becomes the half of it: tr — ^f 2 ^. 
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6 Conclusions and Discussions 



In this paper, we studied a stochastic motion of a uniformly accelerated charged particle in the 
scalar-field analog of QED. The particle's motion fluctuates because of the thermal behavior 
of the uniformly accelerated observer (the Unruh effect). Because of this fluctuating motion, 
Chen and Tajima [TJ conjectured that there is additional radiation besides the classical Larmor 
radiation. On the other hand, it was argued [U [5] that interferences between the radiation field 
induced by the fluctuating motion and the quantum fluctuation of the vacuum may cancel the 
above additional radiation. The cancellation was shown in the case of an internal detector, but 
it was not yet settled whether the same kind of cancellation occurs in the case of a fluctuating 
charged particle in QED. 

In the present paper, in order to investigate the above issue systematically, we first for- 
mulated a motion of a uniformly accelerated particle in terms of the stochastic (Langevin) 
equation. By using this formalism, we showed that the momenta in the transverse directions 
actually get thermalized so as to satisfy the equipartition relation with the Unruh temperature. 
Then we calculated correlation functions and energy flux from the accelerated particle. Par- 
tial cancellation is actually shown to occur, but some terms still remain. Hence there is still 
a possibility that, besides the classical Larmor radiation, we can detect additional radiation 
associated with the fluctuating motion caused by the Unruh effect. 

There are several issues to be clarified. First in calculating the energy flux at infinity there 
appeared classically unacceptable contributions (i.e. those depend on r+). If the observer is 
in the right wedge, the contribution to the energy flux come from the particle in the future of 
the observer. In the case of the observer in the future wedge, this contribution comes from the 
virtual particle in the left wedge. Both of them are classically unacceptable, and we do not yet 
have physical understanding why these contributions appear in the calculation. 

Another issue is the calculation of longitudinal fluctuations. Since the particle is accelerated 
in the longitudinal direction with very high acceleration, longitudinal fluctuations caused by the 
Unruh effect are technically difficult to evaluate (see Appendix [B]). Furthermore, it is not clear 
what kind of quantities are thermalized in the longitudinal fluctuations. The particle feels finite 
temperature noise in the accelerating frame and we may expect that the kinetic energy written 
in the Rindler coordinate £ is thermalized. However, the particle is uniformly accelerated in 
constant external force. Then if a particle at Pi on a trajectory T\ in Fig. [3] is kicked in the 
longitudinal direction to P2, it follows another trajectory T 2 with a different asymptote. The 
Rindler coordinate £, which is defined in the original accelerating frame, becomes divergent at 
t = 00 on the new trajectory T2. In this sense, the longitudinal fluctuations look kinematically 
unstable in the original Rindler coordinate and it is inappropriate to describe the thermalization 
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of longitudinal fluctuations in the Rindler coordinate £. We give brief discussions on the 
calculations of longitudinal fluctuations in Appendix iBl 




Figure 3: T\ is a trajectory of the original accelerated particle. If the particle is kicked at P\ to 
a point P2, it follows a different trajectory T2 in the constant external force. It has a different 
asymptote from the original one, and the original Rindler coordinate £ of the particle on T2 
diverges at r = 00. 

Finally an interesting possibility is an effect of decoherence induced by interaction with 
environments. In this paper, we have treated the trajectory of the particle semiclassically in 
terms of the stochastic Langevin equation. If the initial state of the particle is a superposition 
of two localized wave-packets, we need to sum over the corresponding trajectories to obtain 
the transition amplitude [19J. Decoherence would suppress the correlations between these 
different worldlines. Furthermore the effect of decoherence might be important even for a 
single trajectory with the stochastic fluctuations. We have discussed the partial cancellation 
in the energy-momentum tensor between the inhomogeneous ( fluctuation originated) terms 
and the interference terms. The fluctuating motions of the particle are quantum mechanically 
induced by the vacuum fluctuations of the radiation field. The inhomogeneous terms in energy- 
momentum tensor are evaluated by calculating the correlation of the vacuum fluctuation at the 
almost same positions on the trajectory, and robust against decoherence. Namely the variances 
of the trajectory never vanish. On the contrary, the interference terms are given by calculating 
the correlation function of the vacuum fluctuations at the position of the observer Oj^i and at 
the particle's position on the trajectory, and can be easily affected by an additional interaction 



26 



in between the trajectory and the observer. Then the (partial) cancellation is lost and only 
the correlation function between the inhomogeneous terms (4>i n h{x)4>inh{y)) (namely, the Unruh 
radiation) may survive. 
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A Internal Detector in 3 + 1 dimensions 

In this appendix, we give a brief review of an internal detector in (3+1) dimensions to see how 
the Unruh radiation is canceled by the interference effects [T5j . 

The action for an accelerated internal detector coupling with a massless scalar field is given 



where d T is used to denote a derivative with respect to the proper time r. The equations of 
motion are given by 



by 




(A.l) 




(A.2) 



(d 2 T + n 2 )Q(r) 



A 



<f>(z(r)). 



(A.3) 



m 



Substituting the solution = <ph + 4>inh, 




(A.4) 
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to the equation of the internal detector, we get the following equation, 

+ njj)Q(T) - - / dT'Q(r')G R (z(r)-z(r')) = -M4r)). (A.5) 
m J m 

Here 0^ is the homogeneous solution representing the vacuum fluctuations. The inhomogeneous 

term is evaluated by expanding the Green function with respect to (r — r') as we did in (12.271) . 

Then after a renormalization of the mass term, we get the diffusive term of the radiation 

reaction, 

dr'Q(r')G R (z(T) - z(t')) 5^. (A.6) 
The stochastic equation can be solved by the Fourier transformation on the path as 

Q(r) = Xh(u)<p{u)), (A.7) 

where 

uj\ 2 

h{u)~ l = -my 2 + mtt 2 - %—— (A.8) 
and the Fourier transformations are defined as 

Q(u) = J dre^ T Q(r), (A.9) 

<p(u>)= [ dre^M^r)). (A.10) 



Note that G r (z(t) — z(t')) is a function of (r — r') if the classical solution z(t) represents the 
accelerated path (12. 8p . The 2-point correlation function is decomposed into 

(00)00)) =((l>h(x)<f> h (y)) + ((j)i n h(x)(j) h (y)) + {(j) h (x)(j) inh {y)) + ((j)inh(x)(j) iTlh (y)) (A.ll) 

where 

(<f)inh(x)(f>h(y)) + (4>h(x)(f>inh(y)) 

dr^-e~^X 2 h(uj){G R (y - z{r)){<j> h {x) V {u)) + G R (x - z{r)){ V {u)^ h {y))) (A.12) 

(4>inh(x)4>inh(y)) 

drdr ,^du^ e _ i{ujT+UJ , T , )x4GR{x _ z{T))GR{y _ z ^)) h ^ u ) h ^)^(u)ip(cy)). (A.13) 



We first evaluate the interference term (IA.12I) ; 



h (x)<p{u)) = / dre^(0 o (x)0o(z(r))) 



dr- 



4vr 2 J (x° - z°(t) - ie) 2 - (x 1 - z 1 ^)) 2 - p 2 

(A.14) 
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Poles are given by solving the equation, 



= I x u - =1 { - U - 2^Y_ p 2 (A.15) 



a 



u Vu hx -. (A.16) 

a a a 2 



The solutions of this equation are classified according to two different types of observers 
(see FigfJ) 

Of (in future wedge) : u > 0, v > 



2u\ V a 



=±(-L 2 + Jl' + ^uv) (A.17) 



-^ t = U- L, -)l Lt + 7 m ) (A ' 18 » 

Or (in right wedge) : u < 0, x° + x 1 > 



e arR = 


a 


2\u\ 


e ar + = 


a 


2\u\ 



L 2 — \l L 4 — -z\uv\ ) (A.19) 



a 2 



(L 2 + Jl*--\uv\), (A.20) 



where, L 2 = —x 2 + 1/a 2 . The poles at t^ ,r correspond to the proper times at the intersections 
of the particle's world line and the past light-cone of the observer's position. Hence they are 
the physically acceptable poles. On the other hand, r+ correspond to the proper time at a 
point on a "virtual path" in the left wedge. t+ lies at an intersection of the world line and the 
future light-cone of the observer. Both of them are classically unacceptable. 
Summing these contributions to the integral, we obtain 

frn 

P M = — ^nri^ -e^Z x (u)), (A.21) 



where 



Using the following relation, 



Z x = e ™ /a 9(u) + 6(-u), (A.22) 
Po =^Jl 4 + ^uv. (A.23) 



drG R (x-z(r))f(r) = J_/( T _) J (A.24) 

47rp 
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a part of the interference term depending on or r_ can be written as 

(M^inhiv)) *>? J drdr'^G R (x - z(r))G R (y - z{r'))e^ '> i . (A.25) 

Similarly, we have 

(<t> inh (x)My)) -> *A 2 J drdr'^G R {x - z{r))G R {y - g (^)) e ^-0 _^L_ ) (A .26) 
where we have used the identity 

(0(u)My)) =((My)<p{-u))Y- (A.27) 

The correlation function of inhomogeneous terms is given by 

=A 4 1 drdr'^^e-^^G R (x - z(r))G R (y - z(r'))h(u)h(-uj) 

x J d{T a -r b )2ii5{u + uy^- T *\MATa))M<n))) 

^ j drdr^e^G R{ x - Z (r))G R (y - z(r'))^^t^. (A.28) 

These three contributions ( 1A.25I) . ( 1A.26I ). ( 1A.28I) to the correlation function are shown to be 
canceled each other because of the relation 

h(u>) -h(-cu) = l -^-\h{u)\ 2 . (A.29) 
2ir 

Therefore if we neglect the contributions from the classically unacceptable poles at r + the 
2-point function vanishes, and therefore there are no energy-momentum flux after the thermal- 
ization occurs. 

The remaining term in the 2-point function is the contributions of the r + dependent terms 
to the interference term, and written as 

— ~ m l X \ . L—^( h f u ) e -™ (r-W-r+M) z _ h{-u)e-^ T +^- T -^ Z x {-u)). 

27r87rp Q (x)po(y)l-e- 2 ™/ aK v ; yK ' v ; xK 11 

(A.30) 

It looks strange why we have such a (classically unacceptable) term in the final result. 
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B Longitudinal fluctuations 



4 = —■ (B.2) 



In this appendix we briefly study fluctuations along the direction of the particle's classical 
motion. It is convenient to define the light-cone coordinates 

z ± = z t ±z x . (B.l) 

The classical solution in the light-cone coordinates is given by 

±e ±ar 

a 

In order to describe the fluctuations around the classical solution, we define the Rindler coor- 
dinates (f, £) by 

p ±a(r±0 

z ± = ±- . (B.3) 

a 

The classical solution corresponds to f(r) = r and £(r) = 0. 

Small fluctuations around the classical solution (IB. 21) are written as 

6z ± = (6T±6£)e ±aT . (B.4) 

Writing the velocities as 8^ = ±v ± e ±aT , the stochastic equations for the longitudinal fluctua- 
tions become 

m?) ± = (ii ± ± av* - aV) + (av T T v T ) + eQOcj) (B.5) 

127T 127T 

where 

O = a + e aT d+ -e~ aT d_. (B.6) 

Here 0(t) acts on 4>(z(t)) and d± = |(9 ± d x ). 

The set of equations can be solved by using the Fourier transformations 

w ± (r) = J ^e-^v±(u), (B.7) 

00 = J ^e-^O^co) (B.8) 

as 

12-7T 

^(w) = — 2 -—eQO V (u). (B.9) 

Since (w + — v~) is not affected by the quantum field </?, we can safely put it at zero. It is 
consistent with the gauge condition z ■ z — 1, i.e. i c ; • 8z — 0. 
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By using the relation 

v ± = 5£ + a5f ± (5f + a5£) (B.10) 

we can obtain fluctuations for f and £ as 

8t(u) = ^r— 2T - 2 -—-eQO^u) (B.ll) 

5f(u) = 2 ° 2 12? " 9 . s eQO<p{u). (B.12) 

cr + cj(e 2 a; — 12z7rm) 

The vacuum noise fluctuations for 0<p can be calculated as 

W(r)O0(r')) = ,) — - ■ (B.13) 



It is interesting that this is exactly the same as the noise correlation (13.71) appearing in the 
stochastic equation for the transverse momenta. However the 2-point function of 5^ and Sf 
behave differently from the 2-point function of the transverse momenta. After symmetrization, 
we have 

(Si(r)5i(r')) 3 =6e 2 Q 2 / du coth (™) (^—X 4 ? + (B.14) 



(5f(T)5f(r'))s =6e 2 Q 2 / coth — , - ' -. (B.15) 

x v ; v ;/ ^ 7 V a / \a 2 + u 2 J e 4 w 2 + (12vrm) 2 v 7 

The structure of the integral is quite different from that appeared in the transverse fluctuations. 
The poles at are the same, but the poles at in the transverse case are replaced 

by poles at ±ia. Since a 3> the longitudinal fluctuations are affected by higher frequency 
modes of the quantum vacuum fluctuations. Hence,the positions of the poles are of the same 
order as a and we cannot use the derivative expansion with respect to u/a, which was used 
in the case of the transverse fluctuations. Because of this, we do not know yet the validity of 
the integrals and an appropriate way to evaluate them. This technical problem will be related 
to another problem mentioned in the discussions that it is not clear what kind of quantities 
are appropriate to describe the thermalization in the longitudinal direction. We leave further 
analysis for future investigations. 
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